Local suppression of the hidden order phase by impurities in URu2Si2 
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We consider the effects of impurities on the enigmatic hidden order (HO) state of the heavy- 
fermion material URu2Si2. In particular, we focus on local effects of Rh impurities as a tool to 
probe the suppression of the HO state. To study local properties we introduce a lattice free energy, 
where the time invariant HO order parameter ip and local antiferromagnetic (AFM) order parameter 
M are competing orders. Near each Rh atom the HO order parameter is suppressed, creating a hole 
in which local AFM order emerges as a result of competition. These local holes are created in the 
fabric of the HO state like in a Swiss cheese and "filled" with droplets of AFM order. We compare 
our analysis with recent NMR results on U(Rha;Rui-i,)2Si2 and find good agreement with the data. 

PACS numbers: 71.27+a, 75.40.Mg, 76.60.-k, 74.62.Dh 



I. INTRODUCTION 



The physics of heavy-fermion materials is fascinating 
and extremely challenging due to a variety of exotic phe- 
nomena that can be observed, e.g., the Kondo effect, 
the heavy mass renormalization, the onset of novel mag- 
netism or of unconventional superconductivity. The in- 
terplay between these phenomena makes a detailed un- 
derstanding of the ground state complicated. Here we 
focus on URu2Si2, the heavy-fermion material that ex- 
hibits a "hidden order" (HO) phase below Tho = 17. 7K^. 
The specific heat of this material displays the typical 
jump of a second order phase transition at Tho, how- 
ever the precise nature of the HO remains a subject of 
intensive debate. Far above the HO transition the mag- 
netic susceptibility has a maximum around T ~ 50K— . 
The measured magnetic moment reported by neutron 
scattering, if there is any, is too small (~ 0.03 fis/U) 
to explain the large entropy loss at Tho within a lo- 
calized AFM scenario, which led to the concept of the 
small moment antiferromagnetism^ . Early /iSR (muon 
spin relaxation) measurements reported magnetic mo- 
ments as small as ^ 10^'^ fis/U^- However, later /iSR 
and nuclear magnetic resonance (NMR) measurements 
on pure URu2Si2 revealed an inhomogeneous coexistence 
between the HO and AFM order with a sizable magnetic 
moment^i^. Moreover recent neutron scattering experi- 
ments evidence that this small moment is not an intrin- 
sic feature of the HO, but a spurious effect due to local 
strains induced by crystal defects in the sample^. Nev- 
ertheless magnetic ordering is not completely extraneous 
to URu2Si2: an antiferromagnetic phase with large mo- 
ment can be stabilized by applying pressure or straii*^!^^— . 
Since 1985 several theories have been proposed to iden- 
tify the nature of the hidden order parameter. Recently 
a resurgence of interest in this material has been seen 
as new data and new ideas on the nature of the HO 
appearediSrJi. 

Further progress in experimental techniques such as 
sample quality and more accurate measurements of 



URu2Si2 suggest that a breakthrough in this long stand- 
ing problem is at hand and may be achieved soon. In this 
paper we focus on the role of impurities as probes of the 
nature of the hidden order puzzle. We address the role 
of deliberately placed Rh impurities on the suppression 
of the HO state. Since few impurities are added to the 
sample, NMR is a particularly useful bulk probe sensitive 
to the local atomic environment to reveal what happens 
to the HO state at the impurity site. Recently the ^^Si 
NMR spectrum has been reported in U(Rui_2;Rha;)2Si2 
as a function of temperature and Rh concentration—. 
The experiment showed local suppression of the HO state 
and the emergence of satellite NMR peaks, indicating 
the onset of local antiferromagnetic droplets near each 
Rh impurity. These experiments were interpreted in a 
Ginzburg-Landau framework, where antiferromagnetism 
and hidden order are coupled through gradient terms. In 
this scenario the antiferromagnetism is not a competing 
order parameter but rather a parasitic effect induced by 
spatial inhomogeneities in the hidden order parameter—. 

Here we turn to a more microscopic description of the 
effects of impurities at the atomic length scale by using 
a lattice free energy, where each lattice site corresponds 
to a uranium atom. We thus extend earlier work, using a 
lattice free energy with parameters describing the phase 
diagram of URu2Si2i^ in presence of pressure and strain, 
in order to address a spatially inhomogeneous setting. 
The approach is not tied however to a specific micro- 
scopic origin of the hidden order: the form of the lattice 
free energy is general and we choose a particular set of 
parameters values since it has been proven to be con- 
sistent with experiment. In this model droplets emerge 
around the impurities as a result of the competition be- 
tween the HO and AFM ordeiVjenhanced by the coupling 
mechanism suggested in Ref. [is!. 

The paper is organized as follows: in Sec. II we intro- 
duce the lattice free energy that we minimize in order 
to determine the phase diagram of U(Rui_a;Rhj;)2Si2; 
in Sec. HI we show our results and compare them with 
recent experimental NMR dataii^; in Sec. IV we draw 
our conclusions. Finally in appendix A we discuss how 
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the lattice free energy can be derived from a microscopic 
Hamiltonian and in appendix B we derive the coarse- 
grained Ginzburg-Landau free energy from our lattice 
free energy model to make contact with earlier worki^. 

II. MODEL 

We write the free energy in terms of two order param- 
eters: the HO parameter ipi and the AFM order parame- 
ter Mi. Mi is the z component of the magnetic moment, 
since it is observed experimentally that URu2Si2 orders 
magnetically along the z direction^i^i^. The free energy 
contains three terms: F — + Fm + Fc, where Fc is 
the coupling term between -0^ and Mi . Assuming that 
the hidden order preserves time-reversal symmetry, the 
simplest form of Fc is Fc = gi ^'f^i ~- Therefore 
we can write the lattice free energy as 

i i ij 

+ aM E Mf + bM E + ^ E J!fM.,M, (1) 

i i ij 

+ g^Y.^jMf , 

i 

where ipi and Mi are defined for each site of a three di- 
mensional lattice. 

The form of this lattice free energy is general and can 
accommodate different scenarios for the hidden order in 
URu2Si2. As discussed in the appendices, the informa- 
tion about the underlying microscopic theory is contained 
in the values of the lattice free energy parameters. A 
similar phenomenological free energy was proposed for a 
toy model describing the competing AFM and hexade- 
capolar order emerging from crystal field splitting within 
the unit cell of URu2Si2^^. In that case, the parameters 
were naturally expressed in terms of an effective crystal 
field splitting A at each uranium site: = aM = cl — 
#coth(A), ^ bM = b = f [sinh(f ) - ^ r:£(A/lT] 
and gi — 2b. The effective exchange constants J'^ 
and J*^ were determined in such a way to reproduce 
the experimentally observed critical temperatures, i.e., 
the hidden order transition temperature at zero doping 
Tho ~ 17.7K and the Neel temperature Tn — 15.7K, that 
is J'f' = A/tanh(A/2Tho) and J*=^ = A/ tanh(A/2rAr)2a. 
This parametrization of the lattice free energy, using the 
measured elastic constants, was shown to provide an ex- 
cellent description of the phase diagram of URu2Si2 un- 
der applied magnetic field, pressure and strain . Here 
A = 35K is the effective crystal field splitting between 
the two lowest lying states of the U atom 5/ electrons 
of URu2Si2 in the paramagnetic phase^^. In this work 
we choose to adopt the same parametrization. We stress 
however that the form of this free energy is much more 
general and describes a situation where ^pi is any order 
parameter that does not break time-reversal symmetry; 



for example alternative order parameters are a charge 
density at incommensurate momenta or a hybridization 
wave as proposed in Refs. [H and[23l. 

To incorporate the role of impurities we consider two 
effects. The first is a mean-field effect in which we re- 
gard the coefficients a and b to be disorder dependent. 
In an itinerant picture the presence of disorder creates 
a random potential acting on the electrons: the impuri- 
ties act as scattering centers which reduce the excitonic 
pairing in the particle hole channel. In the model of 
Ref. [l^ the impurity induced strain increases the crys- 
tal field parameter A and therefore it suppresses both 
the antiferromagnetism and the hidden order stabiliz- 
ing the paramagnetic phase. The coefficients a and b 
thus acquire a linear (at leading order) dependence on 
doping X by imposing A = Aq -I- xAi. However we 
keep the definition of J'^ and J*^ to be disorder in- 
dependent, i.e. J'^ = Ao/tanh(Ao/2Tho) and J^^ = 
Aq/ tanh(Ao/2r;v). With this parametrization there is 




FIG. 1. (Color online) Phase Diagram of the lattice free 
energy described in the text in the , J*^ , T space at doping 
a; = 0. 

a duality between hidden order ijj and large antiferro- 
magnetic moment M. When J"^ > J*^ the hidden order 
phase is stabilized at low temperature; if J^^ > J"^ the 
large moment phase is stabilized. In Fig. [T] we show the 
phase diagram of the homogeneous lattice free energy 
of Eq. [T] with the choice of parameters discussed above, 
in the ( J'^, J*^, T) parameter space. The introduction of 
impurities, through the disorder dependence of the mean- 
field coefficients a and &, suppresses equally both HO and 
AFM order parameters. 

Disorder also has a very local effect: the coupling of the 
U which is closest to a Rh site has a different exchange 
interaction with its neighbors. Close to Rh impurities HO 
is suppressed and U atoms acquire a finite moment. This 
is explicitly described by a locally renormalized exchange 
term in the hole of the swiss cheese fabric 

Fn-^\[j''-j'']Y.^''^^''h.^3 (2) 

d=±l 

where > and by imposing il){ri) =0. As a result, 
away from the impurities J*^ < J'^ and the hidden order 
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is stable, but in the immediate vicinity of the impurity 
antiferromagnetism is stabihzed over hidden order. In 
order to study the local suppression of the hidden order 
by impurities, we consider the lattice model described by 
the free energy 



impurities 



= F 



(3) 



To introduce a minimum number of parameters we 
limit the range of the interaction to only the first neigh- 
bors of the affected uranium site. Given these definitions, 
the free parameters of the model are the magnetic cou- 
pling J^^ at the impurity sites and Ai. Due to the duality 
of the model, as magnetic droplets can be stabilized in 
a hidden order background, with the same mechanism 
droplets of "hidden order" could be stabilized within the 
large moment phase by another type of impurity, which 
would exchange the role of J*^ and J*^. It would be in- 
teresting to see if it is possible to realize this dual scenario 
experimentally. The existence of localized regions of the 
HO phase in the AFM phase at higher pressure could 
be observed in NMR experiments and other local probes. 
The counterpart of the previous URui_a;Rh^Si2 experi- 
ment would require to measure URu2Si2 under pressure 
to stabilize AFM but doped with suitable impurities to 
induce a local expansion in the ab plane. 

In the numerical simulations we will assume that the 
disorder is dilute enough to consider the solution of the 
mean- field equations around a single impurity at tq. In 
this work we take a simplified cubic lattice instead of 
the tetragonal lattice of URu2Si2, and we do not con- 
sider the problem of how the order parameter on the U 
atoms is transferred to the nuclear sites of the Si where 
the NMR is performed. Our goal in this paper is to ex- 
plore the physics introduced by inhomogeneities using a 
lattice free energy framework and see how the NMR ex- 
periments constrain the symmetry and the parameters in 
this theory. We test whether the parametrization of the 
lattice free energy that was used to successfully describe 
the phase diagram of URu2Si2 under pressure, stress and 
applied magnetic field can also account qualitatively for 
the puzzling NMR measurements when impurities are in- 
troduced in the sample. 

Once we define the lattice free energy, we determine 
the value of V'i and Mi that minimize the free energy, i.e. 
the solutions to the equations 



-^impurities 

Hi 
-^impurities 



= 
= 0. 



(4) 



Notice that, away from the impurity, the solutions are 

Mf = and ipf = J'^--iaiT,x) ^^^^ ^^^^ ^ ^^^^ ^^^^^ 



enough distance from the impurity the lattice transla- 
tional invariance is restored). 



III. RESULTS 

We first compute the HO critical temperature as a 
function of doping x. Since we assume dilute doping, 
we find 



Ao -I- xAi 



2artanh(^2i^) 



(5) 



At the critical doping Xc the HO parameter vanishes and 
T\io{xc) — 0; it follows that Xc is given by Aq -I- XcAi — 
. In Fig. [2] we compare the theoretical curve ([5]) 




0.01 0.02 0.03 
X (doping) 



FIG. 2. (Color online) Computed critical temperature Tho for 
the order parameter t/i as a function of doping. Squares and 
circles are experimental values taken respectively from Refill 
and Ref. [24!. The dashed line corresponds to the temperature 
Tafm at which the magnetic signal disappears according to 
our calculation. Inset: computed jump of the specific heat 
Acv/Acv(0) at the phase transition as a function of doping. 
Acv(O) is the value of the jump at zero doping. 



for Ai = 358. 6K with experimental data^®'^^. J"^ and 
Aq have been defined above for the uniform case and 
are equal to J''' = 46.24K and Aq = 35K. We ob- 
serve that the transition temperature, up to leading 
terms, has in our model the usual linear dependence on 
X similar to impurity-averaged theories. In the inset of 
Fig. [2] we report the computed jump in the specific heat 
Acv — Since the free parameter Ai was deter- 

mined by the critical temperature the agreement with 
the experimental data is very satisfactory^. In Fig. [2] we 
plot also Tafm as a function of doping. Tafm is defined 
as the temperature at which the magnetization becomes 
smaller than the minimum observed magnetic moment 

/^0,min 

0.03 ^iB- 

In Fig.[3]we report the NMR frequencies / as a function 
of Rh concentration x. The NMR frequency / is propor- 
tional to the spin moment /Uq at each U site, in particular 
/± = liHo ± A^o)(l + K), where 7 = 8.46MHz/T is 
the gyromagnetic ratio of ^^Si, K = 0.065 is the Knight 
shift. Ho = 7T is the external field, A = 0.36T//xb 
is the hyperfine coupling. In our model we identify the 



4 





0.05 




- 63.1 



0.04 



0.01 0.02 0.03 

X (doping) 



FIG. 3. (Color online) Magnetic moment induced by the 
Rh impurities, in units of /xs (left axis) and in units of the 
corresponding NMR frequency (right axis) as a function of x 
at T = 4 K. The frequencies of the two satellite peaks are 
/± = 'y[H(i ± A/io)(l + K), where 7 = 8.46MHz/T is the 
gyromagnetic ratio of ^^Si, K = 0.065 is the Knight shift, 
Ho = 7T is the external field, A — 0.36T//ifl is the hyper- 
fine coupling and fio is the ordered spin moment of U atoms. 
The red line is the result of our model, while full triangles 
are experimental points taken from Ref. 15. The computed 
magnetic moment is the value of the magnetization at the 
impurity site. 



magnetization at the impurity site Mr^ with /io- The the- 
oretical curve in Fig. El is obtained with j'^' = 33.83K. 
The emergence of a magnetic moment in URu2Si2 doped 
with Rh has been observed also in neutron scattering 
experiments^^. Taking into account that neutrons mea- 
sure the magnetization averaged over volume, the agree- 
ment between the two measures is good. From the NMR 
spectrum we obtain another valuable information: the 
area under the satellite peaks is proportional to the frac- 
tion of antiferromagnetic sites, i.e. the ratio between 
sites with a finite magnetization and the total number 
of sites. The measured antiferromagnetic fraction has a 
non-monotonic behavior as a function of the Rh concen- 
tration x: first it increases linearly with the number of 
impurities, then it reaches a maximum at a: = 0.025 and 
finally it decreases becoming zero after the critical con- 
centration has been reached^. The neutron scattering 
result is consistent with this observation: the intensity 
of the magnetic Bragg peaks has a non-monotonic be- 
havior as a function of x - . Our model offers a simple 
explanation of this non-monotonic behavior. In Fig. 2] 
a real-space representation of the magnetization Mi at 
each site is reported for impurity concentrations x — 0.01 
and X — 0.025. The sites around the impurity develop 
a finite moment. The moment at each site decreases as 
the distance from the impurity increases. The magneti- 
zation is strongly suppressed and ipi recovers the bulk so- 
lution value within few lattice sites, see profile picture in 
Fig. m In the following we will refer to the magnetic sites 
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FIG. 4. (Color online) Upper panels: magnetization density 
in the plane z = at T = 4K for Rh concentration x = 0.01 
(left) and x — 0.025 (right). We display the magnetization in 
the interval [0.002, 0.22]: brighter colors correspond to higher 
values. The red contour corresponds to the magnetic droplet 
boundary as defined in the text. Lower panels: profile of the 
magnetization Mi (blue dashed line) and of the HO order 
parameter ^l^i (black solid line) along the direction (0,y,0) of 
the lattice. Red dots correspond to the intersection of the 
magnetization with droplet boundaries. 



around the impurity with the term "droplet" . We define 
the droplet boundary in such a way that the magnetiza- 
tion of every site inside the droplet is large enough to be 
observed experimentally. We consider that the minimum 
observed magnetic moment is equal to /lo.min — 0.03 fiB- 
We can see in Fig. 2] that the size of the droplet is af- 
fected by disorder. The two competing effects of disor- 
der are evident: on the one hand the number of mag- 
netic droplets increases with the number of impurities, 
on the other hand in our model the size of each droplet 
decreases with increasing disorder. This leads to the 
observed non-monotonicity in the experiments. Notice 
that, in a model where HO and AFM order are cou- 
pled through a term M^jV'f/'P, the behavior of magnetic 
droplets in function of doping is similar^?-. In order to 
put this analysis on more quantitative grounds, we opti- 
mize the free energy for different values of temperature 
and Rh concentration and then we compute the fraction 
of sites with Mi ^ {Mi > //o.min)- Since we made the 
assumption that magnetic droplets are disjoint with av- 
erage spacing d ~ l/x^^^ > 360 nm (here / is the lattice 
constant), we can define the antiferromagnetic fraction 
in the following way 



AFM fraction 



riin X 7T.ii] 



TJin X X , 



(6) 
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where riin is the number of sites inside the droplet, riimp is 
the number of impurities and iVtot is the total number of 
sites. In figure EJ upper panel, we plot the AFM fraction 



disappearance of the magnetic droplet for a given Rh con- 
centration. We observe that at large doping Tafm follows 
the behavior of T^o as a function of x, see Fig. O 
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FIG. 5. (Color online) Upper panel: comparison between 
the measured (green circles) and the computed (red squares) 
AFM fraction, as defined in the text, for different values of 
Rh concentration at temperature T — 4K. The black line is 
a guide to the eye. Lower panel: computed AFM fraction 
as a function of temperature T for different values of the Rh 
concentration x. 



as a function of doping. The curve has the characteristic 
non-monotonic behavior of the experimentally observed 
AFM fraction, which highlights the two competing effects 
of disorder: the chemical pressure and the suppression of 
order. Crucial to this observation is the duality between 
■0 and M, and hence the fact that Rh impurities suppress 
both order parameters. While the sawtooth profile is a 
consequence of our lattice model, there is good agreement 
between the results of the minimization procedure and 
experiment. The sawtooth profile appears because the 
number of sites with Mi > /io,min is a step function of x. 
In fact lowering the temperature below T/v the first site 
to be magnetized is the impurity site rg , then the nearest 
neighboring sites, followed by the next-nearest neighbors 
and so on. Therefore the droplet size increases (and with 
the same mechanism the droplet decreases as a function 
of doping) in steps equal to the coordination number. In 
the lower panel of figure O we plot the AFM fraction as a 
function of temperature. At low temperature the AFM 
fraction is a non-monotonic function of x as discussed 
above; increasing the temperature the magnetic droplet 
can be stabilized only for lower doping values. We iden- 
tify with Tafm the temperature in correspondence to the 



IV. CONCLUSION 

In conclusion, we proposed an analysis within a mean- 
field lattice free energy to reveal the local competition 
between HO and AFM phase. We introduced disorder 
in the model as the driving force of two competing ef- 
fects: the local stabilization of magnetization and the 
suppression of both the HO and AFM order by the im- 
purity. We recovered the main features of the phase dia- 
gram and the non-monotonic behavior of the AFM vol- 
ume observed experimentally. Moreover, we found that 
the healing lengths of tpi and Mi are on the scale of the 
local strains that stabilize the magnetization. An addi- 
tional effect present in the calculation is the stabilization 
of the phases due to inhomogencities pointed out in Ref. 
[Tsl . which is manifest in the continuum theory supported 
by our lattice model as shown in Appendix B. Finally 
our model describes a duality between HO and AFM or- 
der: as magnetic droplets can be stabilized by impurities 
in the HO phase, with a similar mechanism HO droplets 
can be formed in a magnetic phase. 

We used in this paper a classical lattice free energy. 
Since the phase transitions occur at rather low temper- 
atures, it would be interesting to extend our work to in- 
clude quantum effects, including effects which would be 
described by time derivatives of the order parameter as 
well as the effects of damping due to particle-hole excita- 
tions. The study of these effects as well as their deriva- 
tion from microscopic models will provide further con- 
straints on the possible origin of the hidden order state 
in URu2Si2. This would allow also a more refined mod- 
eling of the NMR line-shapes and the mechanism for the 
transfer of the hyperfine fields from the uranium to the 
ligand site. 
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Appendix A: Mapping from microscopic 
Hamiltonian onto the lattice free energy 

A general and formal expression for the lattice free 
energy in terms of the HO and AFM order parameters 
ipi and Mi is given by 
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dAi / dAo / X>^X>^'te-/da;£[*^*]+E, Ai,.(Oi([*+*],)-V'.)+i:.A2,.(02([*t*,])-M,) 

J J (Al) 



-PF[ipi,Mi] 



The Lagrangian is the starting point of the cal- 

culation written in terms of creation and annihilation 
operators containing all the relevant bands. The start- 
ing point can be simplified depending on the itinerant or 
localized model one considers and on the full quantum 
many-body Hamiltonian. However, in the general case 
the evaluation of the free energy is more involved and 
will not be attempted here. ip and M are two order pa- 
rameters with ijj time reversal invariant and M breaking 
time reversal symmetry, the index i runs over the lattice 
positions of the U atoms. Different proposals of the hid- 
den order differ in the definition of the operator Oi . In 
the proposal of Ref. [12I the hidden order is a charge den- 
sity wave with incommensurate wave vector Q*, and in 



this case the condensate is related to the Fourier trans- 
form of Oi([*t^f],) 

Oi([vl/tvi/]) = ^fi{k-Q*)^>,,{k),,S„,,, . (A2) 

k,(T(T' 

For a hexadecapolar order as in Ref. [l^ the operator Oi 
is equal to 

Oi ( [vi/t VI/] ) = c ^ ^ F,,, (fc - Q)*^' (fc') 

kk' (J a' 

(A3) 

where C is a normalization constant, Q is the commen- 
surate ordering vector and Fa-a'{k, k') are defined by 



Fnik,k') 

Fnik,k') 
Fu{k,k') 



647r 
3V5 
647r 
3\/5 
647r 
3^5 
647r 



bV7{k^ ~ lkyf{k% - ik'^k'y - k^)k', - {kl + lk.,ky - kl)kSk'y + k',^){~l + bk kll- 

5V7{k^ + ikyf{k^ + ik'^k'y - ky )k'^ - [kl - ik^ky - kl)k^{k'^ - ik'y){-l + bk^) jklk' 
5{kl - ik^ky - kl)KXk^ - ik'^k'y ~ k';f)k', - (fc, - iky)^{k'^ - zfc;)(-l + 5k^ 7fe7*, 
5{kl + ik,ky ~ kl)k,{k^ + ik'^k'y ~ k^)k', - (fc, + ikyf{k'^ + ik'y){-l + bk^)] -fkYk' , 



(A4) 



with 7fc = 47r / dr j^{kr)R{r). In the definition of 7^, 
R{r) is the radial wave function of the /-electrons and 
j^{kr) is the spherical Bessel function of order 3. The 
operator O2 is the magnetization operator 



02([*^*] 

x{Qx 

y.-^a-{k') 



-2MB 



E 



kk' .a,a' 



dfe 



t'rj 



(A5) 



where a are the Pauli matrices and Q is a reciprocal 
lattice vector. 

From a lattice free energy perspective different micro- 
scopic models result in different values of the coefficients. 
An important coefficient is the coherence length; from the 
numerical simulation we estimate the coherence length 
for the magnetic droplet to be ~ 3 lattice constants at 
T = 4K and doping x = 0.01. The NMR data place 
important constraints on this parameter given that more 
itinerant models give rise to longer coherence lengths and 



more diffuse domain walls for the order parameter defined 
on the lattice scale. 



Appendix B: Mapping from lattice free energies to 
coarse-grained Ginzburg Landau free energy 



The free energy of Eq. ([T]) defines a model on the lattice 
describing the 5/-U electrons in URu2Si2. In the contin- 
uum a Ginzburg-Landau (GL) free energy can be derived 
from a lattice model by suitable coarse graining. The 
GL description is formulated in terms of slowly varying 
amplitude fields (f>{x). Here we describe the coarse grain- 
ing starting from the microscopic free energy F used in 
Ref. [1^. We keep higher order terms in the coupling be- 
tween the hidden order parameter and the magnetization 
but focus only on the form of the GL action to connect 
it to the earlier work of Refs. JJi andi21i. 

The starting point is the free energy ([T]) on a cubic 
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lattice introduced in Ref. [l^ 
F[V^„ M„ hf, hf] = 1 ^ - ^ /if 

ij i 

(Bl)' 

written in terms of the order parameters ipi, Mi and of 



the molecular Weiss fields hf , hf-, see Ref. [IQj. At the 



extrema of the free energy and ipi satisfy the rela- 
tions hf = Y^- jfjipj and tjji = — ■^tanh(Ai/TAi) with 

A, = ^(A/2)2 + (/if)2 + (/if )2 i9. The same equations 

are satisfied by hf'^ and M^. In Eq. (|Bip we write /if^ and 

/if in terms of Mi and "0^ exploiting the above relations 
and expand the free energy neglecting terms of the order 

of O ((E, 4^j? + (E, J!fM,ry and higher. The 
Fourier transforms of the lattice variables tpi and Mi are 
V'(fc) = and M(^) = -^^.e^*^-«^Mi. 

Since our mean-field description includes only nearest 
neighbors antiferromagnetic coupling, the modes that 
condense are M{k = Q) or ip{k = Q) with Q = (tt, tt, tt), 
depending on the relative size of J'^ and J*^. The free 
energy of Eq. (jBl[) can be rewritten in terms of ipi^), 
M{k) and the coupling constants J^^^'^\k). In order to 
obtain the free energy in the GL form we keep only the 
modes with k close to Q^^ . For sake of simplicity we shift 
the wave vectors of the Brillouin zone by Q and there- 
fore we consider only the modes close to fc = (fc < A). 
For small k values, the Fourier transform of the coupling 
constant J"^^^") (fc) can be approximated as 



(B2) 



ber. After writing the free energy in terms of ip{k) and 
M{k) and keeping only the modes with /c < A, we go back 
to a real-space representation using the transformation 



k<A 
^ ^ k<A 



M{k) . 



(B3) 
(B4) 



where we scaled J as J/z, z being the coordination num- 



Indeed in Eqs. (|B3|) and (|B4|) the sum is on a discrete non- 
periodic set of wave vectors k with fc < A, therefore the 
fields 01 (r) and (f>2{^ are slowly varying and continuous. 
In terms of the fields 0i (r) and 02 (r) the free GL energy 
becomes 

F= fdfY. (\^^c.iT)[M^y + \h |v0„(f) 

a=l,2 ^ 
Q,/3 = l,2 

-I k^ Yl \yK{r)\\Mr)f 

a, 13=1,2 
a, 13=1,2 ij 

The coefficients ki and u are definite positive. To ob- 
tain the traditional form of the free energy we restrict 
the temperature dependence to the coefficient fi(T) of 
the quadratic term. The coupling coefficients are the 
same for 0i(r) and 02 (^) since this GL free energy has 
been derived from a microscopic model where hidden or- 
der and magnetization are related to each other, however 
the form of the free energy is completely general, and for 
the URu2Si2 system was first discussed in the work of 
Ref. [21I . For a sufficiently repulsive quartic interaction, 
it captures the competition and interplay between the 
HO and AFM order: the field 02(^0 can develop only 
if the hidden order field 0i(r) is suppressed. The gra- 
dient coupling term /c2(02(?^)^|V0i(^)P was introduced 
and discussed in detail in Ref. |l8| to explain the non- 
monotonic behavior of the antiferromagnetic fraction in 
the NMR spectrum. When k2 < inhomogeneities in 0i 
can nucleate a parasitic second order parameter 02 near 
impurities even when /^2 > 0. 



^ T. T. M. Palstra, A. A. Menovsky, J. vandenBerg, A. J. 

Dirkmaat, P. H. Kes, G. J. Nieuwenhuys, and J. A. My- 

dosh, Phys. Rev. Lett. 55, 2727 (1985). 
^ C. Broholm, J. K. Kjems, W. J. L. Buyers, P. T. Matthews, 

T. T. M. Palstra, A. A. Menovsky, and J. A. Mydosh, 

Phys. Rev. Lett 58, 1467 (1987). 
3 D. E. MacLaughlin, D. W. Cooke, R. H. Heffner, R. L. 

Hutson, M. W. McElfresh, M. E. Schillaci, H. D. Rempp, 

J. L. Smith, J. O. Willis, E. Zirngiebl, C. Boekema, R. L. 

Lichti, and J. Oostens, Phys. Rev. B 37, 3153 (1988). 



^ G. Luke, A. Keren, L. P. Le, Y. J. Uemura, W. D. Wu, 

et al, Hyperfine Interact. 85, 397 (1994). 
^ K. Matsuda, Y. Kohori, T. Kohara, K. Kuwahara, and 

H. Amitsuka, Phys. Rev. Lett 87, 087203 (2001). 
^ P. G. Niklowitz, C. Pfleiderer, T. Keller, M. Vojta, Y. K. 

Huang, and J. A. Mydosh, Phys. Rev. Lett. 104, 106406 

(2010). 

^ M. Yokoyama, H. Amitsuka, K. Tenya, K. Watanabe, 
S. Kawarazaki, H. Yoshizawa, and J. A. Mydosh, Phys. 
Rev. B 72, 214419 (2005). 



8 



A. Villaume, F. Bourdarot, E. Hassinger, S. Raymond, 
V. Taufour, D. Aoki, and J. Flouquet, Phys. Rev. B 78, 
012504 (2008). 

^ E. Hassinger, G. Knebel, K. Izawa, P. Lejay, B. Salce, and 

J. Flouquet, Phys. Rev. B 77, 115117 (2008). 
^° C. R. Wiebe, J. A. Janik, G. J. MacDougall, G. M. Luke, 

J. D. Garrett, H. D. Zhou, Y.-J.Jo, L. Balicas, Y. Qiu, 

J. R. D. Copley, Z. Yamani, and W. J. L. Buyers, Nat. 

Phys. 3, 96 (2007). 
" A. F. Santander-Syro, M. Klein, F. L. Boariu, A. Nuber, 

P. Lejay, and F. Reinert, Nat. Phys. 5, 637 (2009). 

A. V. Balatsky, A. Chantis, H. P. Dahal, D. Parker, and 

J. X. Zhu, Phys. Rev. B 79, 214413 (2009) 
" P. Chandra and P. Coleman, Nature Physics 5, 625 (2009). 

F. Cricchio, F. Bultmark, O. Granas, and L. Nordstrom, 

Phys. Rev. Lett 103, 107202 (2009). 

K. Haule and G. Kothar, Nature Physics 5, 796 (2009). 

S. Elgazzar, J. Rusz, A. M. Amft, P. M. Oppeneer, and 

J. Mydosh, Nature Materials 8, 337 (2009). 

C. M. Varma and L. Zhu, Phys. Rev. Lett. 96, 036405 

(2006). 

S.-H. Baek, M. J. Graf, A. V. Balatsky, E. D. 
Bauer, J. C. Cooley, J. L. Smith, and N. J. Curro, 



Phys. Rev. B 81, 132404 (2010) 

K. Haule and G. Kotliar, EuroPhysics Letters 89, 57006 
(2010). 

'° C. Broholm, H. Lin, P. T. Matthews, T. E. Mason, W. J. L. 
Buyers, M. F. Collins, A. A. Menovsky, J. A. Mydosh, and 
J. K. Kjems, Phys. Rev. B 43, 12809 (1991). 
N. Shah, P. Chandra, P. Coleman, and J. A. Mydosh, 
Phys. Rev. B 61, 564 (2000). 

Notice that the coefficients and J*^ contribute both 
to the massive and gradient term in the continuum for- 
mulation. In fact in order to derive the continuous free 
energy, we substitute Jij = J{R) with the Fourier trans- 
form J{k) = EHe''°'^J(-R) (see Appendix B). The zero 
order term of J{k) contributes to the massive term, while 
the quadratic term in k to the gradient term. 
2^ P. Woffle, Y. Dubi, and A. V. Balatsky, Phys. Rev. Lett. 
105, 246401; Y. Dubi and A.V. Balatsky, unpublished. 
(2010). 

M. Yokoyama, H. Amitsuka, S. Itoh, L Kawasaki, 
K. Tenya, and H. Yoshizawa, J. Phys. Soc. Japan 73, 
545 (2004). 

D. J. Amit, Field Theory, the Renormalization Group, and 
Critical Phenomena (World Scientific, 1984). 



